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ABSTRACT 

Extreme ultraviolet (EUV) lithography is seen as a main candidate for production of future generation computer 
technology. Due to the short wavelength of EUV light (w 13 nm) novel reflective masks have to be used in the 
production process. A prerequisite to meet the high quality requirements for these EUV masks is a simple and 
accurate method for absorber pattern profile characterization. 

In our previous work we demonstrated that the Finite Element Method (FEM) is very well suited for the simu- 
lation of EUV scatterometry and can be used to reconstruct EUV mask profiles from experimental scatterometric 
data. 13 ' 16 

In this contribution we apply an indirect metrology method to periodic EUV line masks with different critical 
dimensions (140 nm and 540 nm) over a large range of duty cycles (1:2, ... , 1:20). We quantitatively compare 
the reconstructed absorber pattern parameters to values obtained from direct AFM and CD-SEM measurements. 
We analyze the reliability of the reconstruction for the given experimental data. For the CD of the absorber 
lines, the comparison shows agreement of the order of lnm. 

Furthermore we discuss special numerical techniques like domain decomposition algorithms and high order 
finite elements and their importance for fast and accurate solution of the inverse problem. 

Keywords: EUV, mask, metrology, photolithography, FEM 

1. INTRODUCTION 

The quality of pattern profiles of EUV masks becomes important e.g. due to shadowing effects at oblique 
incidence illumination. 1 ' 2 Consequently, there is a need for adequate destruction free pattern profile metrology 
techniques, allowing characterization of mask features down to a typical size of 100 nm and below. The desired 
accuracy is thereby of the order of 1 nm which is challenging even for direct atomic force (AFM) or scanning 
electron microscopic (SEM) measurements. These methods e.g. suffer from surface charges and the need of the 
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Figure 1. (a) Experimental setup for scatterometry experiment with fixed incident angle of 8i n = 6° and variable angle 
of detection 8 ut', (b) Result of single wavelength scatterometry measurement at A = 13.655 nm. Diffraction orders appear 
as peaks with finite width, the zeroth diffraction peak is centered around 6°. 

definition of edge operators to reconstruct the CD. Furthermore it is extremely time consuming and not practical 
to scan large areas of a mask with direct microscopical measurements. 

Here the usage of scatterometry as a metrology technique offers the possibility to analyze larger areas on a 
mask depending on the diameter of the incident beam, e.g. arrays of periodic absorber lines, contact holes etc. 
Of course detecting small errors of a single feature remains extremely challenging. 

2. CHARACTERIZATION OF EUV MASKS BY EUV SCATTEROMETRY 

Single wavelength scatterometry, the analysis of light diffracted from a periodic structure, is a well suited tool 
for analysis of the geometry of EUV masks. Since scatterometry only needs a light source and a simple detector 
with no imaging lens system, its setup is inexpensive and offers no additional technical challenges. Fig. 1(a) 
shows a sketch of the experimental setup. Light of fixed wavelength and fixed incident angle is reflected from 
the mask and the intensity of the reflected light is measured in dependence on the diffraction angle. The use 
of EUV light for mask characterization is advantageous because it fits the small feature sizes on EUV masks. 
Diffraction phenomena are minimized, and of course the appropriate wavelength of the resonant structure of the 
underlying multilayer is chosen. Light is not only reflected at the top interface of the mask but all layers in the 
stack contribute to reflection. Therefore one expects that EUV radiation provides much more information on 
relevant EUV mask features than conventional long wavelength methods. The presented EUV measurements 
often provide up to 20 or more non-evanescent diffraction orders. 

All measurements for the present work were performed by the Physikalisch-Technische Bundesanstalt (PTB) 
at the electron storage ring BESSY II. 3 PTB's EUV reflectometer installed at the soft X-ray radiometry beamline 
allows high-accuracy measurements of very large samples with diameters up to 550 mm. ' 5 

Experimental data was taken from an EUV test mask fabricated by AMTC. The mask is subdivided into 11 
(labeled A-K) times 11 (labeled 1-11) test fields. Each of these test fields furthermore consists of a number of 
differently patterned areas allowing e.g. bright field, dark field and CD uniformity measurements. In the present 
work fields with periodic absorber lines of varying pitch and duty cycle are considered. We focus on fields D4, 
H4, F6, D8, H8 (uniformity fields with 720 nm pitch and 540 nm CD) and fields HI to H5 (varying pitch and 
140 nm CD), see Table 1. 

3. FEM SIMULATION OF EUV SCATTEROMETRY 

Fig. 1(b) shows the result of a scatterometry measurement of an EUV test mask (parameters given in Table 1) 
considered in the present work. 
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Table 1. Design parameters (see also Fig. 4) for EUV test mask and test fields produced by AMTC. 




Figure 2. SEM pictures of EUV mask patterns and corresponding triangulated geometries for FEM computation. 

The positions of the diffraction angles provide information about the pitch of the EUV absorber pattern. 
However the intensities of the diffraction orders do not carry direct information about other topological features 
of the mask. The determination of these features from a scatterometry measurement is a challenging inverse 
problem and a hot topic of actual research. Accurate and fast numerical simulation of the scattering experiment 
thereby plays a vital role. The experimental situation of single wavelength scatterometry is described by time- 
harmonic Maxwell's equations for fixed frequency oj. We choose the finite element method to obtain solutions 
to Maxwell's equations numerically. The idea of the finite element method is very simple and will shortly be 
explained in the following. Maxwell's equations can be transformed without approximations into weak form. 
Therefore Maxwell's equations for the electric field E are multiplied by a so called test function 3? and integrated 
over the computational domain fi (e.g. the EUV mask): 

/ {$ • [V x x E - oj 2 e ■ E] } d 3 r = 0, (1) 

Jn 

where the bar denotes complex conjugation. With a partial integration (Stoke's theorem) this equation can be 
transformed to 

J |(V x *) • (/i" 1 V x E) - w 2 e ¥ • E j d 3 r = J * ■ Fd 2 r, (2) 

where F contains data of the electric field on the boundary T = dfl. Defining the complex bilinear form a and 
linear form /: 

a(#,E) = J {(V x *) • (/^V x E) - w 2 eI .E}d 3 r, (3) 
/(*) = j^-¥d 2 r, (4) 



Maxwell's equation can be stated compactly in weak formulation (i.e. as a variational problem) 6 : 



Find E £ V = H(curl, fl) such that: 



o(#,E) = /(#), V*GV. (5) 

H(curl, f2) is an inhnite dimensional function space and the exact solution to Maxwell's equations is an element 
of this space. The exact mathematical definition of H(curl,Q) and the theory behind the weak formulation 
can be found e.g. in. 6 The finite element method discretizes the infinite dimensional space V by restricting 
it to a finite dimensional subspace Vh C V . The weak formulation (5) is then stated on this finite dimensional 
subspace and therewith discretized. For the performance of the finite element method for Maxwell's equations 
it is very important to construct Vh in such a way that the geometrical properties of V carry over to Vh- The 
appropriate construction with so called edge elements e.g. assures the continuity of the solution at boundaries of 
different materials, a fundamental property of the electric field (in the usual case of absent surface charges). The 
finite clement method does not approximate Maxwell's equations itself (a(., .) and /(.) are not changed) only the 
the solution space is approximated with a finite dimensional version Vh- The construction of Vh is done in the 
following. First the computational domain is subdivided into small patches, e.g. triangles in 2D or tetrahedrons 
in 3D. On these patches a local finite elements space is defined which usually consists of polynomial functions (so 
called shape functions) up to a certain order p. Since the electric field is vector valued the shape functions are 
also vectorial. To increase the accuracy of a finite element solution either the size of the patches can be reduced 
[h refinement) or the polynomial order can be increased (p refinement). For wave propagation phenomena as 
described by Maxwell's equations "high order" is often a significantly better strategy then "fine grid" as we will 
see in a convergence analysis at the end of this section. 

It has been shown that the FEM method is an excellent choice for this type of application. 7 It produces 
very accurate results in relatively short time which is vital for the solution of the inverse problem. The FEM 
method has several advantages 7 : 

• Maxwell's equations describing the scattering problem are solved rigorously without approximations. 

• The flexibility of triangulations allows modeling of virtually arbitrary structures, as illustrated in Fig. 2. 

• Adaptive mesh-rcfincment strategies lead to very accurate results and small computational times which 
are crucial points for application of a numerical method to the characterization of EUV masks. 

• Choosing appropriate localized ansatz functions for the solution of Maxwell's equations physical properties 
of the electric field like discontinuities or singularities can be modeled very accurately and don't give rise 
to numerical problems. 

• It is mathematically proven that the FEM approach converges with a fixed convergence rate towards the 
exact solution of Maxwell-type problems for decreasing mesh width of the triangulation. Therefore it is 
easy to check if numerical results can be trusted. 

Throughout this paper we use the FEM package JCMsuite developed by the Zuse Institute Berlin and JCMwave 
for numerical solution of Maxwell's equations. JCMsuite has been successfully applied to a wide range of 
electromagnetic field computations including waveguide structures, 8 DUV phase masks, 7 and other nano- 
structured materials. 9 ' 10 It provides higher order edge elements, multigrid methods, a-posteriori error control, 
adaptive mesh refinement, etc. Furthermore a special domain decomposition algorithm implemented in JCMsuite 
is utilized for simulation of EUV masks. 11 Light propagation in the multilayer stack beneath the absorber pattern 
can be determined analytically. The domain decomposition algorithm combines the analytical solution of the 
multilayer stack with the FEM solution of the absorber pattern, decreasing computational time and increasing 
accuracy of simulation results. The iterative idea of the domain decomposition algorithm is depicted in Fig. 4(b). 
In the first step light scattering from the absorber pattern is computed with the finite clement method neglecting 
the multilayer (1). The resulting field on the lower boundary is decomposed into diffraction modes (2) and the 
reflection of each diffraction mode from the multilayer is computed analytically (3). In the next step the incident 
field from above and the reflected diffraction modes from below are coupled into the computational domain 
(i.e. the absorber pattern) and the FEM computation is repeated, yielding a new solution for the electric field. 
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Figure 3. Relative error of complex diffraction orders A(/ n ) in dependence of number of unknowns of FEM computation 
for different diffraction orders n, finite element degrees p and refinement strategies, (a) Regular refinement, FEM degree 
p = 6; (b) regular refinement, diffraction order n = 15; (c) finite element degree p — 6. 



This solution is again decomposed into diffraction modes on the lower boundary and the iteration is continued. 
This iterative solution converges to the exact solution of the coupled problem, 11 i.e. domain decomposition 
as implemented in JCMsuite is a "rigorous" method. 12 For the problem class investigated in this paper, few 
iterations are enough to reach a relative error < 0.1% for the diffraction orders. It is worth noting that the FEM 
system has to be assembled only once. Also the LU-factorization (i.e. inversion) of the FEM matrix is only 
performed once and applied to the new right hand side in each iteration step. 



4. CONVERGENCE OF FEM SIMULATIONS 

Fig. 3 shows the convergence of the complex diffraction orders for the simulated EUV mask with Layout H5, 
see Table 1. The error is given relative to the most accurately computed solution. Since we look at the complex 
diffraction orders phase errors are also included. Fig. 3(a) shows the relative error of the diffraction orders 
n = —9, 0, 5, 15 in dependence on the number of unknowns for regular refinement and finite element degree 
p = 6. All orders are converging at almost the same rate. Fig. 3(b) shows the big benefit of higher order 
elements. For first order finite elements (degree p = 1) we see almost no convergence, which is in agreement with 
the slow convergence of first order finite difference time domain (FDTD) methods for photo mask simulations,. 7 
For orders p = 4 and p = 6 the convergence rate of the diffraction orders increases dramatically. The flexibility 
of the FEM method allows to choose optimal settings of grid refinement h and polynomial order p of the ansatz 
functions to achieve a specified accuracy in minimum computational time. Fig. 3(c) shows a comparison of the 
convergence rate for uniform and adaptive refinement for finite clement degree p = 6 and the diffraction orders 
n = 10, 15. Adaptive refinement increases the convergence rate even more. For the same number of unknowns 
the relative error is almost 2 orders of magnitude smaller. In the following computations we used finite clement 
order p = 4, and a grid refinement level where the relative error of the diffraction orders is below 0.1%. 



5. INVERSE SCATTEROMETRIC PROBLEM 

Determining the profile of the periodic absorber lines we have to solve the following optimization problem: 

Optimization problem: For given experimental diffraction orders I exp = (In Xp )\ n= i n determine the geom- 
etry G of the EUV mask such that: 

G = arg mind{I exp ,r imG '), (6) 



are the diffraction orders obtained from the scattering simulation with 

n=l,...,JV 



where 1*™*°' = (l» im >°'^ 

geometry G' . The metric d measures the difference between experimental and simulated diffraction orders. 
S is the set off all possible geometries. 




Figure 4. (a) Geometrical parameters describing absorber line profile, (b) Illustration of scattering process from EUV 
mask: (I) incident plane wave, (2) light is decomposed into diffraction orders by absorber line, (3) each diffraction order 
is reflected with amplitude and phase shift according to multilayer properties, (4) scattered light from EUV mask. Of 
course also light coming from the multilayer is reflected again by the absorber pattern and also part of the incident light 
is reflected directly from the absorber pattern. 



As stated above the optimization problem is infinite-dimensional and ill-stated. The set of possible geometries 
S is infinite-dimensional and we only have a finite number N of diffraction orders used for determination of G. 
One can imagine that there arc a lot of different geometries which produce similar diffraction patterns. 



The following questions and problems arise solving (6): 



• Experimental diffraction orders I exp are measured with uncertainty. 

• The material properties (permittivity) are not known exactly. 

• The layout of the multilayer beneath the absorber line is not known exactly. 

• How do we restrict 2 to arrive at a better conditioned problem? 

• Which error measure do we use for d ? 

Having these points in mind we can not expect that we find a geometry G rea i such that d(I exp , I si ™< G ™i ) = 
and probably a lot of geometries produce errors d of similar magnitude. Consequently the question of accuracy 
and error bounds for the minimizing geometry G op t of (6) arises. 

The experimental diffraction orders are measured with a relative accuracy of 1%. Noise in the EUV detector 
is of the order of 5 ■ 10 -5 relative to the incident light. Experimental diffraction orders with intensity below this 
limit are therefore not used for optimization since their relative error is above 100%. For d we use the sum of the 
squares of the relative errors over all diffraction orders. In the present work experimental orders were taken at 
three different wavelengths for each absorber pattern. All of the orders (with intensity > 0.5 • 10 -5 ) are included 
computing d: 



For definition of the set of geometries S considered in minimizing (6) we parameterize the absorber line as 
depicted in Fig. 4(a). We have three discrete parameters: the height of the absorber stack h, the absorber 
sidewall angle a and the width of the absorber line, e.g. given at the top, i.e. top CD. Here we will only look at 
the results for the CD and a. The pitch p of the periodic lines can be deduced from the position of the diffraction 
peaks. The subset of possible geometries S/, C 3 is therewith defined. 

Accurate modeling of the multilayer is a crucial part in the reconstruction of the absorber pattern. 13 We think 
that in our setting it has the greatest influence on reconstruction results and is the main cause for uncertainty 
of the reconstructed parameters. Other error sources are e.g. numerical errors of the FEM simulations. They 
can be virtually neglected using appropriate computational settings after a convergence analysis, according to 
Section 4. In 14 the uncertainty of the reconstruction results caused by measurement errors of the diffraction 
orders is analyzed and quantified, which is also important to judge the reliability of reconstructed geometry. 

The experimental orders of our analysis have a measurement uncertainty of 1%. The diffraction orders of the 
numerically determined best fitting geometries however still differ much more from the experimental orders then 
these measurement errors, see Fig 9. Consequently at this stage we focus on the reconstruction error caused by 
inaccurate modeling of the multilayer. 



The scatterometric measurements of the EUV mask were performed close to the later operating wavelength of 
the EUV mask. Fig. 5(a) shows a typical reflectivity curve of a blank multilayer stack. A maximum reflectivity 
of 60% to 70% is reached when the component of the fc-vector of the incident plane wave which is orthogonal to 
the layers is in resonance with the half period of the multilayer. These so called bright field curves were measured 
on unstructured regions close to each of the test fields of the EUV mask. Each scatterometric measurement was 
performed for the left and right full-width half-maximum (FWHM)-wavelength and center wavelength X C enter 
(mean of right and left FWHM wavelength), see Fig. 5(a). 

The usage of resonant EUV light has a great advantage. The intensity of reflected experimental diffraction 
modes is increased significantly since the EUV stack reflects these modes well in the resonant case. This makes 
it possible to measure more diffraction orders above the limiting intensity of the EUV detector. Especially the 
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6. MODELING OF THE MULTILAYER 
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Figure 5. (a) Different numerical fits of experimental bright field curve of field H4. For each curve a different wavelength 
interval of the measured bright field curve was fitted: left to center wavelength, around center wavelength and right 
to center wavelength, (b) Bright field measurement and diffraction orders n of field H4 in dependence on wavelength. 
Intensity of diffraction orders are scaled for better comparison. The wavelength A = 13.655 nm corresponds to the center 
wavelength of the bright field curve. 



high diffraction orders carry the most information about the absorber pattern 13 and are therefore important for 
a well conditioned reconstruction. 

Light scattering from an EUV mask can be decomposed according to Fig. 4(b). It follows that the intensities 
of the diffraction modes used for reconstruction depend on the shape of the absorber pattern and the multilayer, 
which is demonstrated in Fig. 5(b). The wavelength dependency of the diffraction orders resemble the bright 
field curve to a great extend. The disadvantage of using resonant EUV light is obvious. If the multilayer is not 
modeled correctly the absorber pattern can not be reconstructed accurately. Therefore experimental data (the 
bright field curve) is needed and has to be fitted by modeling the multilayer before reconstruction. 

Since the multilayer curves for all test fields had a similar shape only bright field curves from two different 
fields (H4 and H5) were used modeling the multilayer. For both of these curves the multilayer profile was 
computed by fitting simulated reflectivity to the experimental values. As material stack a Molybdenum/Silicon 
(MoSi) multilayer with MoSi2 intcrlayers was used. Diffusive or rough interfaces between the materials were not 
taken into account. Optimizing a complete measured reflectivity curve over a large wavelength interval was a 
difficult task and the results were not very satisfactory for the EUV test mask. Difficulties in the modeling are 
the mentioned interlayers and the unknown thickness profile of the multilayer. When growing a multilayer it is 
hard to determine how well the nominal values of the multilayer are reached. Thinking of a multilayer with e.g. 
50 bilayers with inter- or diffusion layers of varying thickness one arrives at a optimization problem with several 
hundreds of unknowns. Consequently modeling of the multilayer itself is a difficult task and actual research 
problem. 15 

Here we decided to model three small wavelength intervals of the complete multilayer bright field curve 
independently, left to the center wavelength, around the center wavelength, right to center wavelength. In total 
this gives 6 different multilayer curves (3 for H4 and H5 respectively) which were used for reconstruction of the 
absorber profile. The fits of the bright field curve for field H4 are shown in Fig. 5(a). The "fit: right" curve fits 
the right side wiggles of the experimental curve best, the "fit: left" the left side wiggles. Although fitting only a 
small wavelength interval we still see rather poor agreement for the side wiggles for all fits. 

In total we get 6 reconstructed geometries, one for each modeled multilayer. These are compared to quantify 
the uncertainties of the reconstruction process. Although the experimental bright field curves of the considered 
test fields had a similar shape their center wavelengths differed slightly, see Table 2. Scaling all layers of a 



modeled multilayer equally, each of the 6 bright field curves used for reconstruction could be shifted to match 
the center wavelength of the field it is used for. We also compare reconstruction results with and without this 
shift. 

Test field center wavelength [nm] Test field center wavelength [nm] 



HI 13.582 D4 13.6479 

H2 13.615 H4 13.6510 

H3 13.639 F6 13.6769 

H4 13.655 D8 13.6470 

H5 13.664 H8 13.6518 



Table 2. Center wavelength of bright field curves of test fields. 



7. RECONSTRUCTION RESULTS 

We applied the described reconstruction procedure to the test fields with nominal values given in Section 2. The 
profiles of the absorber lines were measured with direct AFM and CD-SEM for a quantitative evaluation of the 
indirect reconstruction results. The microscopically measured values are given in table 3 and 4. We notice that 
the CD of the uniformity field of H4 differs significantly from that of the other fields although the nominal value 
of 180 nm is the same for all fields. Fig. 6 shows the explanation, the lines for the uniformity field H4 are not 
etched to the bottom. We can already expect poor reconstruction results for this field. 16 

test field CD [nm] CD uniform, [nm] test field CD [nm] CD uniform, [nm] 



HI 138.3 2.5 D4 186.7 3.2 

H2 134.3 4.0 H4 168.0 7.3 

H3 127.8 3.8 F6 180.3 2.7 

H4 129.8 7.1 D8 188.7 3.0 

H5 130.4 5.5 H8 184.1 3.6 



Table 3. SEM measurements of line profiles. CD uniformity is defined as 3<r over an ensemble of measurements, the CD 
is the mean over the ensemble 



test field CD [nm] (thr. A) CD [nm] (thr. B) a[°] (rising) a 2 [°] (falling) 



HI 169.0 160.8 78.7 83.8 

H2 160.9 157.5 87.2 86.5 

H3 154.7 152.2 87.2 86.9 

H4 155.9 153.7 87.6 86.9 

H5 155.5 153.8 87.5 85.9 



Table 4. AFM measurements of line profiles. Sidewall angles were measured in rising and falling mode. The CD was 
measured at two different AFM thresholds. For uniformity fields no AFM measurements were performed due to the small 
bright line width. 



First we look at the reconstruction results of the uniformity fields (D4, H4, F6, D8, H8) with 180 nm bright 
and 540 nm dark lines. Fig. 7(a) shows the range of sidewall angles we obtain when using the 6 different 
multilayer curves with scaling to match X ce nter of each test field (see Section 6). The sidewall angles differ over 
a large range. Since the bright lines of the uniformity fields were to small for AFM measurements we have only 
experimental results for the absorber angle of fields HI to H5, Table 4. However it is reasonable to assume that 
the sidewall angle of the uniformity fields is similar. With this assumption the reconstructed sidewall angles of 
the uniformity fields differ significantly from the experimental values. Fig. 7(b) shows that also the top and 
bottom CD differ significantly from the CD SEM values and except from test field F6 and D8 using 6 different 
multilayers produces CDs over a large range. Field H4 shows the greatest variations as expected because the 



Figure 6. Top down SEM image of uniformity fields D4 and H4. 
completely etched. 

(a) (b) 



It can be seen that the lines for field H4 are not 




Figure 7. Minimum and maximum (a) sidewall angles and (b) top and bottom CDs obtained from reconstruction with 
6 different multilayer curves. 



bright lines were not etched through to the bottom, Fig. 6. At first glance the results look very disappointing. 
With a sidewall angle ^ 90° the CD of an absorber line however depends on the height in which we measure it 
CD = CD(h). The question is if we can find a height h op t where the reconstructed CDs are in best agreement 
with the microscopically determined SEM CDs. Fig. 8(a) depicts the mean CD error p (the mean was taken 
over the 5 test fields) between reconstructed and SEM CDs in dependence on the height h for the 6 modeled 
multilayers. The mean error p is defined as: 

P(h) = #tcst 1 fields E I ( CD i EM ~ (CD SEM )) - (CD( n "(h) - (CD^(h))) | , (8) 
test field i 

where CDf EM is the CD determined by SEM and CD{ nv {h) the CD at height h of the absorber line de- 
termined by inverse reconstruction. The mean value (...) is taken over all test fields of the summation, i.e. 
{HI, H2, H3, HA, H5} and {DA, HA, F6, D8, H8} respectively In the definition of p a constant offset between 
both methods is neglected. Also the CD measured by AFM and SEM differ by an offset. Looking at Fig. 8(a) we 
see that although producing wrong top and bottom CDs the reconstructed CDs at the height of 40 nm to 50 nm 



agree very well with the SEM CDs for all used multilayers. This is slightly above the half of the total absorber 
height. For each multilayer h opt (which minimizes p) was determined. It would be desirable to analyze more test 
fields, e.g. with absorber lines of different heights to check if h op t is always at the same position h opt relative 
to the total height of the line. If that would be the case one could use this value h op t for all reconstructions 
and would not need prior CD SEM measurements to minimize p and therewith determine h op t- We take the 
mean of the six obtained h opt values: h opt = (h opt ). Of course an analysis of a greater number of test fields 
could give a better estimation for h opt . For a quantification of the uncertainty of the reconstructed CDs we now 
determine the CDs at height h opt = (h opt ) from our simulations with 6 different multilayers. To clarify again: for 
each of the 6 modeled multilayers we get a different h opt , which can only be determined using the microscopical 
results for the CD. Since in practice we want to use the indicrect method without SEM measurements we have 
to choose h op t- Our experience has shown to choose it slighlty above the half absorber height. With fixed h opt 
the CDs from the simulations with the different multilayers will not be in best agreement with the SEM results 
and we will utilize this disagreement for our error estimations. The minimum, maximum and mean CDs of the 
6 simulations with different multilayers at same fixed h op t are shown in Fig. 8. The reconstructed results and 
error bounds (minimum and maximum CD at h opt ) are well within the error bounds of the CD SEM results 
which demonstrates the robustness of the reconstruction of the CD very nicely. Field H4 again shows larger 
disagreement. Surprisingly the values obtained from the indirect reconstruction even for this field seem more 
reasonable than the CD SEM results, Fig. 6 shows that the lines are not etched through but still have same CD. 

To summarize the sidewall angles of the absorber lines can not be determined very well due to inaccurate 
modeling of the multilayer, however the CD at the middle of the absorber line can be reconstructed very well. 
Fig. 9 shows the diffraction orders which were computed for the best fitting geometries for two example test 
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Figure 8. (a) Mean error p (8) of CD (over all uniformity test fields) between inverse reconstruction and SEM measure- 
ments in dependence on absorber height h, from which the reconstructed CD was taken, (b) Minimum maximum and 
mean CD at (h op t) from indirect reconstruction with 6 different multilayers, CD SEM values with CD uniformity (with 
error bars). 



fields. Since we used 6 different multilayers we also get 6 intensities for each diffraction mode. Fig. 9 shows the 
minimum and the maximum of these diffraction intensities. We notice that almost all experimental diffraction 
orders are located between the smallest and largest simulated diffraction orders. A better modeling of the 
multilayer is therefore the next step for improving the reconstruction performance, e.g. for the sidewall angle. 
The measurement uncertainties are still much smaller than the tube spanned by the simulated orders, and 
therewith it is confirmed that the measurement uncertainties are of minor importance. 

For the 140 nm lines we compared the results obtained with and without scaling of the multilayer with 
respect to the central wavelength. The results for the CD at h opt are shown in Table 5. Here we performed 
the reconstruction with only one fixed multilayer model. With scaling the relative errors are below 0.5%. The 
absolute error is below 0.7 nm. To our knowledge this is the most accurate CD reconstruction for periodic EUV 
line masks. Without scaling the relative error of the reconstructed CDs is large and > 1% for some fields. 




Figure 9. Experimental and simulated diffraction orders of reconstructed geometry. The minimum and maximum 
diffraction order over all 6 modeled multilayers is given; (a) field H4 with incident wavelength = left FWHM wavelength, 
(b) field F6 with incident wavelength = center wavelength. 

test field CD SEM CD 5 ""' (Agp), with scaling CD 5cgtf - (A C p), without scaling 

~H1 lS3 136.5 (1.3%) 138.9 (0.4%) 

H2 134.3 133.1 (0.9%) 134.0 (0.2%) 

H3 127.8 129.4 (1.3%) 127.9 (0.1%) 

H4 129.8 131.1 (1.0%) 129.3 (0.4%) 

H5 130.4 130.4 (0.0%) 130.4 (0.0%) 

Table 5. Comparison of CD-SEM measurements and reconstructed CDs (given in nm) with scaled and unsealed multilayer 
at height h op t (height with maximal correlation of CD to SEM CD). The relative error Acd is given in brackets. 



8. CONCLUSIONS 

We have shown that EUV scatterometry in combination with FEM simulations is a very well suited method for 
profile reconstruction of EUV masks. A comparison of CD values obtained by inverse reconstruction and direct 
microscopical SEM measurements shows an agreement below 1 nm. The diffraction intensities of a EUV mask 
depend on the actual absorber structure and the properties of the underlying multilayer. Without an accurate 
modeling of the multilayer not all parameters of the absorber profile can be reconstructed, e.g. the sidewall 
angle. However the performance of the reconstruction of the CD is very good and robust. By modifying the 
multilayer uncertainty estimations in the reconstruction results of the CD can be obtained. Those are comparable 
to uncertainties of the CD SEM measurements. 

Differences in the diffraction intensities between the best fitting geometries and the experimental values can 
also be explained by the difficulties of modeling the multilayer with high accuracy. This will be the next step in 
improving reconstruction results. 
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